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Using the measured optical response and surface roughness topography as inputs, we perform 
realistic calculations of the combined effect of Casimir and electrostatic forces on the actuation 
dynamics of micro-electromechanical systems (MEMS). In contrast with the expectations, roughness 
can influence MEMS dynamics even at distances between bodies significantly larger than the root- 
mean-square roughness. This effect is associated with statistically rare high asperities that can be 
locally close to the point of contact. It is found that, even though surface roughness appears to have 
a detrimental effect on the availability of stable equilibria, it ensures that those equilibria can be 
reached more easily than in the case of flat surfaces. Hence our findings play a principal role for the 
stability of microdevices such as vibration sensors, switches, and other related MEM architectures 
operating at distances below 100 nm. 

PACS numbers: 03.70.+k, 68.37.Ps, 85.85.+J 



I. INTRODUCTION 

Electromagnetic fluctuations that pervade any medium 
including empty space generate forces between neu- 
tral bodies known as Casimir-Lifshitz forces, of which 
van der Waals forces are special cases Casimir 
forces arise from electromagnetic waves created by 
quantum and thermal fluctuations These are ex- 
pected to become important as components of micro- 
electro-mechanical systems (MEMS) enter submicrome- 
ter separations! 11 ' 23 ~ — The small scales at which MEM 
engineering is now conducted have revitalized interest in 
the Casimir force since devices such as vibration sensors 
and switches are made with parts that are just a few mi- 
crometers in size and have the right size for the Casimir 
force to play a role: they have surface areas big enough 
but gaps small enough for the force to draw components 
together and lock them tight, which is an effect called 
stiction. In fact, permanent adhesion is a common cause 
of malfunction in MEMS devices. Casimir forces, in syn- 
ergy with electrostatic actuating forces, can further aug- 
ment this phenomenon^ - — Additionally, as the develop- 
ment of MEMS evolves toward nanomechanical systems 
(NEMS), attention will also be drawn to scaling issues. It 
is inevitable that the Casimir interactions between metal- 
lic and/or dielectric surfaces in nanometer proximity of 
each other will occur, and stiction phenomena require 
specific attention. On the other hand, the irreversible ad- 
hesion of moving parts resulting from Casimir and elec- 
trostatic forces j^ 3 - - — can also be exploited to add new 
functionalities to micromcchanical architectures. 

Casimir forces that are responsible for stiction in dry 



conditions^ and thus profoundly influence the actuation 
dynamics, supplement the electrostatic force in counter- 
ing the elastic restoring force to determine, for example, 
the beam's actuation behavior in microswitches. The 
latter is typically constructed from two conducting elec- 
trodes of which one is hxed and the other is suspended 
by a mechanical spring governed by Hooke's lavs4£ (see 
Fig. Hall . Voltage application between the electrodes 
moves the electrodes toward each other because of the 
electrostatic force. At a certain voltage, the moving elec- 
trode becomes unstable and collapses or pulls-in onto the 
ground electrodej 31 i 32 Residual stress and fringing field 
effects have also shown to have a great influence on the 
behavior of RF (radio frequency) switches, and strongly 
influence their failure characteristics ! 33 : 34 

In earlier investigations of the effect of the Casimir 
or van der Waals forces on the dynamical behavior of 
nanoscale electrostatic actuators, roughness was either 
ignored or only weak roughness was considered. In some 
cases tabulated optical data were taken into account^ - — 
In this paper, we will explore the actuation dynamics 
of microswitches made from real materials (with a def- 
inite measured optical response 41 - - — and characterized 
by some degree of nanoscale roughness) accounting for 
both electrostatic and Casimir forces, which counter- 
act an elastic restoring force (see Fig. [Taj. Advances 
made in the measurement and theoretical understand- 
ing of Casmir forces over the last 10 years allow today a 
more detailed study of MEMS made from real material 
surfaces^ - — Note that although electrostatic forces can 
be switched off if no potential is applied, Casimir forces 
will always be present and will influence the actuation 
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dynamics. 



II. ROUGHNESS CORRECTION TO THE 
ELECTROSTATIC FORCE 

It has been show n 44 i 45 that the disagreement between 
the experiment^ and the theory describing roughness 
perturbativel y 48 ' 49 can be resolved by taking into account 
rare high peaks on rough surfaces. These high peaks 
can be described by "extreme value statistics" as follows 
from a statistical analysis of AFM topography data for 
gold filmsj 44 ' 45 i 50 Indeed, recently there has been more 
awareness of the importance of extreme value statistics 
for the analysis of rough surfaces^ The Casimir force 
between rough surfaces can be written as: 



Fcas(z) = F PT (z) + F peaks (z). 



(1) 



The term Fpt(z) denotes the Casimir force between 
rough surfaces from Ref. [48], which includes a "zeroth" 
order contribution corresponding to flat surfaces and a 
perturbative roughness correction ~ (w/z) 2 , where w 
is the root-mean-square roughness. The term F pea ). s (z) 
represents the contribution of high peaks, which is asso- 
ciated with the aforementioned extreme value statistics. 
It is important to note that Fc a s (z) is singular at the dis- 
tance upon contact 44 i 45 i 50 (z = do) which is the real min- 
imum separation due to surface roughness. It is assumed 
that the contributions of high peaks are independent of 
each other, an approximation justified by the large hor- 
izontal distance between them. This distance is large 
because such peaks are statistically rare events i 44 ' 45 The 
ellipsometry measurements for gold samples reported in 
Ref. Hl| were used as optical data. 

The roughness correction to the electrostatic force can 
be obtained in the same way as was done in Ref. |3 for 
the Casimir force: the heights of the surface comparable 
to w can be taken into account perturbatively, whereas 
the contribution of high peaks can be approximated by 
treating each peak independently. This approximation 
is justified by the large distance between the high peaks, 
because such peaks are statistically rare events. The per- 
turbative roughness correction to the electrostatic force 
is based on an analysis for isotropic roughness ; 25 ' 26 ' 46 as 
is the case for the gold films considered here^i which are 
grown under non equilibrium conditions. This correc- 
tion starts by modeling the surfaces as a capacitor with 
capacitance: 



(C(z)) = 



2(2tt) 4 



J(\h(k)\ 
o x 



(k 2 + 



f coth(/tz))dk| (2) 



where the first and second terms correspond to flat sur- 
faces and a second order perturbative correction, respec- 
tively, and A denotes the surface area of each plate. The 
quantity k c represents the wavenumbcr corresponding to 



a lower lateral roughness cutoff of the order of the inter- 
atomic distances (~ 4 A for gold) . For power law or self 
affine random roughness a suitable model for the power 

spectrum ^|/i(k)| 2 ^ to perform calculations with is given 

byS 



|Mk)| 5 



w 2 e 



(2tt) 5 {l + ak 2 e) 1+H ' 



(3) 



Here £ denotes the correlation length, a represents the 
self- affine roughness parameter which can be found by 
solving the algebraic equation: a = l/(2iT)[l — (1 + 
ak 2 ^ 2 )~ H ]. For the gold films considered here^i the 
roughness exponent has the value H = 0.9. 

The electrostatic force including the second order per- 
turbative roughness correction can be written as 



F pe {z) = -\V' 



d{C(z)) 



(4) 



where the average capacitance (C(z)) is given by Eq. @ 
and V denotes the applied voltage. Equipotential planes 
are expected to be a valid approximation at separations 
below 100 nm. Statistical deviations from this approxi- 
mation, known as potential patches, typically play a role 
at separations of the order of a few hundred nanometers 
up to a few micrometers £2r£!i Now the contribution of 
the high peaks in the surface can be approximated by 
a sum of separate contributions of each peak as it was 
done in Ref. for the Casimir force. For this purpose 
we start with the electrostatic force between flat surfaces: 



Fe{z) 



epAV 2 
2z 2 



(5) 



For the roughness statistics we can use the same AFM 
topography data with the same statistical analysis as in 
Ref. [4J. Therefore the electrostatic force between rough 
surfaces becomes 

F es (z) = F pe {z) + (6) 
Jf(d)[F e (z - d) - F e (z) + dF' e {z) - ±d 2 F>'(z)}dd 



where f(d) denotes the probability density function. The 
height d\ = 3w is the separation above which f(d) can be 
fitted to a Gumbel distributio n 44 ' 45 and do is the height 
of the highest asperity (see Fig. [Taj) . It must be noted 
that the expression in Eq. © is also singular at z = do. 

Figure [lb] shows the relative strength of electrostatic 
and Casimir forces for various potentials between real 
nanoscalc rough Au-Au surfaces. The Casimir force be- 
comes significant for separations z < 100 nm and over- 
comes the electrostatic force rather rapidly as the applied 
potential drops below 1 V (a regime typical for MEMS) 
and separations close to distance upon contact due to 
surface roughness. Indeed, the potential V eq (z) where 
Fcas(z) = F es (z) increases rapidly toward smaller sepa- 
rations, which shows that the Casimir force corresponds 



3 



to increasing values of the applied voltage V. These re- 
sults clearly show that below 100 nm Casimir forces can 
strongly influence the actuation dynamics. 

III. ACTUATION DYNAMICS OF MEMS 

Modeling a MEMS as a classical mass-spring system 
has been well established— Let the separation z depend 
on time and satisfy the following differential equation: 

d 2 z 

m-^ = k{L q - z)- F{z), (7) 

where F(z) = Fc< a s(z) + F es (z) represents the total sur- 
face force, k is the spring constant, and Lq is the distance 
between bodies if no external force is present, F(z) = 0. 
The effective mass m merely rescales the equation ([7]). 
For our calculations we used as an example a resonance 
frequency u>o = y n/m = 300 • 2ir rad/s which is typical 
for a wide variety of commercial resonators, e.g. tapping 
mode AFM cantilevers and other doubly clamped beam 
MEMS^. This frequency is kept constant, whereas k is 
varied and used as a control parameter. 

The solutions of Eq. (0 can be investigated with a 
phase portrait^ i.e. a plot of z versus z'(t). Studies of 
the influence of the Casimir force for nanoscale electro- 
static actuators with flat, perfectly conducting electrodes 
showed that their phase portraits exhibit periodic orbits 
around a center equilibrium, and an orbit that passes 
through an unstable saddle point i 38 i 39 These studies were 
extended to the influence of weak roughness onlyi 25 i 26 

Consider first the case of zero electrostatic force, 
F es (z) = 0. The goal is to find out under what condi- 
tions the oscillator described by Eq. (JTJ can return to its 
original position; i.e. for what parameter values periodic 
solutions exist. The existence of periodic solutions indi- 
cates that the spring is strong enough to prevent stiction. 
If the spring constant is large enough, the stable cen- 
ter around which periodic solutions exist will be accom- 
panied by an unstable saddle point equilibrium— - — If 
the spring constant becomes lower, the center and saddle 
point will merge into an unstable "center-saddle" point. 
For an even lower value of k there are no equilibria at 
all. This is an example of what is known as a saddle- 
node bifurcation— 

In order to understand for what values of k such equi- 
librium points are available we introduce the following 
bifurcation parameter 

Acas = L r — , (8) 

where Fl(Lq) denotes the Casimir force given by the Lif- 
shitz formula 2 -! 3 - (for flat surfaces) at z = Lq. This ratio 
of the minimal Casimir force and the maximal elastic 
restoring force represents the relative importance of one 
force compared to the other. In an equilibrium the to- 
tal force given by Eq. ([7]), is zero: F to t{z) — 0. This 



case yields A cas = (1 - z s /L )F L (L )/F cas (z s ), where z s 
denotes the locus of the stationary points. 

The results are plotted in Fig. [2a] As one can see 
the rougher the sample (i.e. the higher the value of the 
contact distance^ do), the higher the spring constant 
must be to get equilibria and periodic solutions. The 
maximum of A cas decreases with do. The position of 
this maximum changes only slightly under the influence 
of random roughness: from 0.78Lo for a flat surface to 
0.81-Lo for the roughest sample. This is because at these 
separations the roughness effect is small (perturbativc) 
and does not drastically change the force. To clarify the 
meaning of Fig. [2aJ the general solution, represented by 
the phase portrait, is plotted for three different values of 
the spring constant, for the roughest sample (with do = 
50.8 nm). Fig. I2bl shows the case where the spring con- 
stant is large enough for the bifurcation parameter to 
be below its maximum value. In this case there are two 
equilibria: the stationary point closest to Lo is a (stable) 
center around which periodic solutions (closed curves) 
exist. Since the system considered here is conservative 
(Hamiltonian), the phase portraits can be obtained by 
plotting the level curves of the total energy. The so- 
lutions of Eq. ([7]) are periodic if the amplitude stays 
below a value of approximately 0.4io- The shift of the 
minimum separation due to roughness from zero to do 
prevents periodic motion if the total energy is too high. 
In Fig. [2c] the value of k has been chosen such that it 
corresponds to the maximum value of A cas in Fig. I2al 
In this case there is only one equilibrium, known as a 
center-saddle point^ which are always unstable. There 
are no periodic solutions in this case. If the value of the 
spring constant is lowered further, no equilibria are avail- 
able anymore. The spring is too weak to counterbalance 
the attractive Casimir force. The solution for this case is 
plotted in Fig. [2d] 

Although neglecting the electrostatic force can provide 
some insight, this force must also be taken into account. 
If we consider the presence of the electrostatic force only, 
F es (z) ^ and Fc as (z) = 0, we can define an additional 
bifurcation parameter for the electrostatic force: 

e AV 2 FejLp) 

Aes = o t3 ~ T • W 

Similarly to the previous case, this is the ratio of the 
minimum electrostatic force and the maximum elastic 
force, which is a measure of the relative importance of 
one force compared to the other. In this case it holds that 
A es = F e (L>o) I F es {z s )(\ — z s /Lq), which is obtained from 
the condition F to t = 0. However, it must be stressed that 
this case is a rather artificial one, because the Casimir 
force cannot be shut down (since it stems from quantum 
mechanical uncertainty). Results for this case are qual- 
itatively similar to the previous one, but the roughness 
effect is less pronounced because the electrostatic force 
depends more weakly on the separation distance than 
the Casimir force. 

In the more general case F es (z) ^ and Fc as {z) ^ 0, 
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the stationary points z s satisfy the following equation 
obtained from the condition F to t = 0: 

-. z s Fcas{ z s) , Fes( z s) , _ „ ,.. „\ 

X -T - YdU) Xcas ~ fm Aes - °' (10) 

where A cas and A es arc defined by Eqs. © and © respec- 
tively. Eq. (|10[) is an implicit function of two variables, 
plotted in Fig. [3a] for both the idealized case of flat sur- 
faces and the roughest surface with do = 50.8 nm. The 
graph for the case of flat surfaces encloses the one for the 
rough surface case. This indicates that, similarly to the 
previous case, surface roughness increases the minimum 
value of the spring constant required to compensate for 
both the electrostatic and Casimir forces. However in 
this case this value further increases with the value of 
the applied voltage V, as indicated by the contour plots 
for the rough case in Figs. [3b] and [He] The phase por- 
traits are similar to those in Fig. [5] the only difference is 
that the distance between the center and saddle point is 
smaller than in Fig. [2a] since V ^ in this case. From 
Fig. [3b] it can be concluded that this distance decreases 
with V since X es cc V 2 . 

The critical equilibrium points at which stiction oc- 
curs are characterized by the conditions dU/dz = and 
d 2 U/dz 2 = 0, where U denotes the total potential energy 
of the system (i.e. F tot {z) = —dU/dz). See e.g. Refs. 
l36l . I38I . or [i^. Hence Eq. ([TO]) must be combined with its 
derivative with respect to z s , which is also set to zero: 

L Fl(Lq) F e (L ) Aes ~ U ' (Uj 

which corresponds to d 2 U/dz 2 = —dF tot /dz = 0. Eqs. 
(fTU]) and (jlip form a system of two equations with three 
unknowns. Eliminating z s from this system yields a re- 
lation between the bifurcation parameters X cas and X es . 
This relation is plotted in Fig. I3dl Since A es oc V 2 / n 
and A cas oc 1/k the origin of this graph corresponds to 
high spring constants and low voltages, i.e. to the situa- 
tion that the spring is strong and the device is shut down 
(i.e. V decreases such that A es < 1). In Fig. [3d]it can be 
seen that the equilibria join in the origin both for the flat 
and the rough case. This corresponds to the ideal case 
where only periodic motion exists and stiction is of no 
concern. It turns out that this limit is reached faster as 
a function of the parameter values in the rough case than 
in the flat one. On the other hand, stable and unstable 
equilbria are closer in the rough than in the flat case. 
It should be kept in mind, however, that the rough case 
in Fig. [3d] corresponds to the roughest sample available 
with do = 50-8 nm. Presumably the optimum in terms 
of preventing stiction lies somewhere between these two 
cases (e.g. a contact distance between 20 nm and 30 nm). 

IV. CONCLUSIONS 

It may seem surprising that a flat surface is not optimal 
in terms of preventing stiction, but the reason for this is 



that the shift of the minimum separation from zero to e?o 
prevents both the Casimir and electrostatic forces from 
becoming too large. In other words, there is a trade-off 
between two effects: on one hand random surface rough- 
ness requires a higher spring constant for equilibria to 
exist and puts the stable and unstable equilibria closer 
together, and on the other hand it prevents the surfaces 
from reaching too short a separation, reducing both the 
Casimir and electrostatic forces. 

Using measured material and surface properties and 
realistic device dimensions, we studied the effect of 
the Casimir-Lifshitz force on the actuation dynamics of 
MEMS. It has turned out that this force is equivalent to 
a voltage between 0.2 and 0.6 V at separations below 100 
nm, which is comparable to the electrostatic force typi- 
cally used to actuate MEMS. We have shown that ran- 
dom roughness has an overall strong effect on the avail- 
ability of equilibria of MEMS oscillators: the rougher 
the surface, the lower the maximum values of the bifur- 
cation parameters, and hence the higher the spring con- 
stant must be in order for equilibria to exist. Finally, the 
shift of the minimum separation due to surface roughness 
moves the stable center curve closer to the unstable sad- 
dle point curve in the bifurcation diagram as a function 
of the stationary points. However, there is a trade-off: 
a MEMS of which the actuating components have rough 
surfaces benefits more from a strong spring constant than 
one with a flat surface. This is because surface roughness 
brings the two equilibrium points closer together in terms 
of the separation, but it prevents the surface forces that 
give rise to stiction from becoming too large. Most likely 
the optimum lies somewhere in between a flat surface 
(with do = 0) and a very rough surface (with do 50 
nm). Filially, another effect that deserves further con- 
sideration is that of nonlinearity, i.e. deviations from 
Hooke's la w 59 ' 60 at sufficiently large amplitudes. For 
example, for operation of a resonator in vacuum stud- 
ied here, we could have some manifestation of nonlinear 
behavior for oscillation amplitudes -~100 nm^i which is 
comparable to the device dimensions, i.e. Lq considered 
here. Surface roughness could possibly reduce or prevent 
such nonlinearities, because high surface peaks limit large 
amplitude oscillations. Our analysis applies to motion in 
vacuum or dilute gases, where friction losses can be ig- 
nored. Qualitatively, the center equilibrium in Fig. [2b] 
becomes an inward 'sink' spiral in the presence of friction. 
A more detailed analysis of the effect of hydrodynamic 
drag forces^r— present in air will be combined with that 
of nonlinearity in a future study. 
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FIG. 1: (a) Schematic of the MEM system. The spring tries to move the oscillator toward the initial separation Lq. The 
contact distance do is the maximum height of the asperities on a rough surface. |(b) | In the inset, the Casimir force from Eq. 



([T]) for gold surfaces is compared to the electrostatic force between rough surfaces from Eq. (J6]) for several values of the applied 
voltage V . The quantity V eq indicates to what voltage the Casimir corresponds after equating the Casimir force from Eq. ((TJ) 
to the electrostatic force from Eq. (J6]). In this case the r.m.s. w — 10 nm and the contact distance do = 50.8 nm. 
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rough vs flat surfaces K = 0.5 N/m, g> /2ji = 300000 Hz 



flat surface 




FIG. 2: (a) Bifurcation diagram of Casimir actuated MEMS for various rough surfaces, each of which is associated with a 
different value of the contact distance do- The value of the initial separation Lq — 100 nm. (b) Phase portrait for the surface 
with do = 50.8 nm. The (black) circle and the (black) square indicate the positions of the center and saddle point equilibria, 
respectively. In this case, the spring constant is high enough for periodic solutions to exist, (c) This phase portrait corresponds 
to the maximum of the solid (red) curve in Fig. I2al In this case there is only one (unstable) equilibrium, and there are no 
periodic solutions, (d) Phase portrait corresponding to a point above the maximum of solid (red) curve in Fig. I2al There are 



no equilibria in this case. 
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FIG. 3: (a) 3D Bifurcation diagram for a MEMS actuated by both Casimir and electrostatic forces for both a rough and a flat 
surface. The innermost graph represents the rough case. |(b)| Side view of the innermost graph of Fig. I3al|(c)| T op view of the 
innermost graph of Fig. I3al (d) Two parameter bifurcation diagram under the conditions of Eqs. (|10[) and . In the origin 
the equilibria join. 



